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PACS 03 . 65 . W j - State reconstruction, quantum tomography 
PACS 71 . 70 . Jp - Nuclear states and interactions 

Abstract. - Precisely knowing an interaction Hamiltonian is crucial to realize quantum informa- 
tion tasks, especially to experimentally demonstrate a quantum computer and a quantum memory. 
We propose a scheme to experimentally evaluate the spin-spin interaction for a two-qubit system 
by the weak measurement technique initiated by Yakir Aharonov and his colleagues. Furthermore, 
we numerically confirm our proposed scheme in a specific system of a nitrogen vacancy center in 
diamond. This means that the weak measurement can also be taken as a concrete example of the 
quantum process tomography. 



Introduction. — Since the publication of the Shor 
algorithm [1], that the prime factorization based on the 
current cryptography scheme can be solved with poly- 
nomial gates in quantum computation, many researchers 
have tried to realize a quantum computer in some physi- 
cal systems [2J. The realization of quantum computation 
is equivalent to designing the Hamiltonian adjusting to the 
given problem. Therefore, we have to precisely know a spe- 
cific spin-spin interaction in a two-qubit system to realize 
quantum computation. We emphasize that the specific 
spin-spin interaction is quite different from an effective 
spin-spin interaction derived from the collective spins like 
the mean-field approximation. From the viewpoint of real- 
izing quantum information processing tasks, knowing the 
effective Hamiltonian is futile since we cannot perfectly 
control a considered qubit. 

Also, quantum storage devices, especially a quantum 
memory, is needed to realize quantum information pro- 
cessing tasks [3|. The vital properties of them are to be 
strong against decoherence, to write, and to read out eas- 
ily. To construct them, we basically prepare two differ- 
ent qubit systems. One is easily manipulated and read 
out, and the other is insensitive to decoherence, i.e., its 
quantum information can be stored for a long time. Fur- 
thermore, we need to transfer quantum information for 
one qubit to the other, which is called to a quantum state 
transfer. There is the problem to exactly know the specific 



spin-spin interaction to perfectly realize the quantum state 
transfer. In a typical example of the Hamiltonian estima- 
tion, the dominant term of the Hamiltonian is assumed. 
Therefore, when we estimate the relative position between 
two spins from the measurement result, the Hamiltonian 
can be evaluated, e.g., the ESR technique [I]. Since de- 
coherence and an experimental error are inevitable in all 
physical systems on measuring the relative position be- 
tween spins, it is difficult to perfectly transfer one qubit to 
the other for the considered two-qubit system. Therefore, 
we need an evaluation of the considered spin-spin interac- 
tion without knowing the relative position to improve an 
accuracy of quantum state transfer. 

A quantum memory with small storage has been exper- 
imentally demonstrated in the systems between the po- 
larization of a single photon and a single-electron spin in 
quantum dots |5j and between the polarization of a single 
photon and the A-type atom in the cavity quantum elec- 
trodynamics (cavity QED), which is called an electromag- 
nctically induced transparency (EIT) [5] . The polarization 
of photons is used as the controlled device of quantum 
communication and quantum computation. The electron 
spin in quantum dots and the A-type atom in the cavity 
QED are used as the quantum storage devices. Also, a 
nitrogen vacancy (NV) center in diamond is a utmost po- 
tential example of the quantum storage devices. A NV 
electron spin state can be mapped to a nearby 13 C nu- 
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clear spin state [7J, which has a long coherence time at 
room temperature. Recently, a NV electron spin state has 
been also mapped to another one, which was located near 
the considered NV center, to improve the scalability of 
this system [5]. 

In this Letter, we propose a scheme to experimentally 
evaluate an unknown specific spin-spin interaction in a 
two-qubit system. We numerically confirm our proposal 
using an example of the NV center in diamond. 

Our ultimate goal is precise parameter estimation of the 
specific spin-spin interaction, which is essential to realize 
quantum information processing tasks, especially a quan- 
tum memory, as discussed above. This is similar to the 
motivation of ref. j9] while they discussed this under dif- 
ferent regulations. The key point of our proposal is to use 
the weak measurement [10] introduced as follows by the 
different viewpoint to the original scheme. 

By taking the weak measurement and the post-selection 
for a target system, we can obtain the weak value of an 
observable A defined as (A) w :— (f\A\i)/{f\i), where \i) 
and (/| are called a pre- and post-selected state, respec- 
tively [10]. These concepts are summarized in refs. [TT|n2] . 
This tool is useful to deeply understand quantum founda- 
tions, e.g., to resolve the Hardy paradox [13] and to give 
an outlook on the macroscopic realism using the viola- 
tion of the Leggett-Garg inequality [Uj . Therefore, many 
theoretical proposals to measure the weak value have re- 
cently proposed in some physical systems, e.g., ref. [T5] . 
However, the problem of the weak value and the weak 
measurement is a little practical advantage, e.g., the first 
demonstration of the spin Hall effect of light [16] and 
the feedback control to stabilize the Sagnac interferom- 
eter [T7]. While the interaction between the target and 
probe systems is given in the conventional concept of the 
weak measurement, our addressed problem is that the in- 
teraction Hamiltonian is unknown and has to be evalu- 
ated p. In our proposal, the short-time evolution can be 
guaranteed in the weak measurement regime. In this Let- 
ter, we also show that the weak measurement can also be 
taken as a concrete example of the quantum process to- 
mography. Therefore, our proposal gives a practical out- 
look on the weak measurement and the weak value. 

Proposed Protocol by Weak Measurement. 

Throughout this paper, we assume the followings: 

1. We can prepare a two-qubit system. One spin is called 
a target spin and the other is called a probe spin. The 
total Hamiltonian is given by 



Htot = H t + H p 



Hi, 



(1) 



to simplify the discussion. Here, H t ,H p , and Hi„t 
denote the target spin Hamiltonian, the probe spin 
Hamiltonian, and the interaction Hamiltonian be- 
tween the target and probe spins, respectively. 



1 Our proposed scheme is different from the weak measurement 
tomography; the post-selected state can be evaluated from the pre- 
selected state and the measured weak value I18| . 



2. We can know the single spin dynamics: H t and H p . 

3. We can manipulate and detect the state of the target 
and probe spins. Note that, the direct measurement 
for spins is not necessary. After transferring the quan- 
tum state to another spin, measuring the transferred 
spin can be taken as the detection of the spin state 
as in ref. [7j. 

The above assumptions are realizable for some physical 
systems in recent quantum information technology. 
Let us consider an unknown interaction Hamiltonian: 



Hi, 



J2 9^{<®°p), 

fj.,ve{x,y,z} 



gxx 9xy 9xz 

[9fJ.v\ = | 9yx 9yy 9yz 

9zx 9zy 9zz 



=: (n x n y n z ) . 



(2) 



where erf and a p are the Pauli matrices (p, v — x,y, z) 
for target and probe spins, respectively. The above equa- 
tion represents the spin-spin interaction between target 
and probe spins. Here, the unknown parameters g^ are 
assumed to be a symmetric tensor and have six degrees of 
freedom. In the following, we propose a scheme to evaluate 
the parameters g^ of eq. ([2]) : 

Step 1: We prepare the pre-selected state for the target 
and probe spins denoted as $i = pt (8> p P , where 
Pt = (I + ri ■ <T t )/2 and p p = (I + p • <r p )/2 in the 
Bloch sphere representation. Here, / expresses the 
identity operator. o~ t and o~ p are the Pauli vectors. 
The controllable vectors, r-j and p, are the Bloch vec- 
tor or the Stokes vector for the initial target spin and 
the initial probe spin, respectively. 

Step 2: We wait for the short time 5t. This proce- 
dure can be taken as the weak interaction regime 
by the unknown spin-spin interaction. The state 
<f>2 = e~ iHtotSt &ie iHtotSt is weakly entangled between 
target and probe spins. 

Step 3: We decide the post-selected state for the target 
spin pt = Tr p $2 by undertaking the quantum state 
tomography. It should be noted that the quantum 
state tomography is not completed in a single event. 
We can repeat this procedure since the state $2 is 
unchanged on remaining the conditions by Step 2. In 
principle, we perfectly undertake the quantum state 
tomography for the target spin. 

Step 4: We measure the probe spin in the measurement 
direction q to obtain the expectation value as 



Ex(q-<7 p )=Tr(P(q)p p ), 



(3) 



where p p := Tr t $2, P (q) := q • o p is the projection 
operator, and q is defined in the following step. 
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Hint 


Parameter Sets 


Ti 


P 


q 


5t[[is] 


/ 5 -6.3 -2.9 \ 

-6.3 4.2 -2.3 
\ -2.9 -2.3 8.2 / 

(in MHz) 


(0,0,1) 


(0,0.59,0.81) 


(-0.16,0,0.99) 


0.091 


(-0.48,0.59,0.65) 


(0,0,1) 


(-0.25,0.59,-0.77) 


0.086 


(-0.81,0.59,0) 


(-0.65,0.59,-0.48) 


(0.25,0.59,-0.77) 


0.073 


(0,0,1) 


(0,0,1) 


(-0.99,0,-0.16) 


0.069 


(0.81,0,-0.59) 


(-0.10,0.95,0.29) 


(0,0.81,-0.59) 


0.066 


(0.31,0.95,0) 


(-0.18,0.95,-0.25) 


(0,0.81,0.59) 


0.051 



Table 1: Example of our proposed scheme. We apply our proposal to the NV center in diamond. The target and probe spins are 
the NV electron and nearby 13 C nuclear spins, respectively. Also, the hyperfine interaction is given by eq. (15) in Supporting 
material in ref. [21] but is taken as the unknown interaction to be evaluated. The set of the controlled parameters is decided 
from Fig. [2] It should be noted that the transposition is omitted in the columns of i"i,p, and q due to the limited page while 
they should be described as the row vector. 



Step 5: We remove the contributions from the single spin 
Hamiltonians H t and H p by calculating the reverse 
operation for the interaction time St to the post- 
selected state Yf for the target spin and the modified 
measurement direction q for the probe spin. Because 
of the short interaction time St, the parameters Yf and 
q can be defined as ^HtSt^-iHtSt _. (/ + r/ • a t )/2 



state tomography 



and 



„iH v St 



P(q)e- 



iH p St 



P(q) by the Trotter-Suzuki 



decomposition 19!. 



Step 6: Changing the controllable parameters; the initial 
target spin state r^, the initial probe spin state p, the 
interaction time St, and the measurement direction q, 
we repeat the procedures from Step 1 to Step 5. 

Step 7: Since Steps 2 and 3 can be taken as the weak 
measurement regime, we obtain the expectation value 
in Step 4 under the first order approximation as 



Ex(q ■ a p ) rs q • p 



+2 St 



2 St 



{(qx rv)-p}( r » + r /)- e P 

1+Ti-Tf 



[n M • q - (n M • p)(q • p)](r, x 17) • e p 



1 



r t 



(4) 



where e M is a unit vector in the y, direction (fi = 
x,y,z). Note that, the above equation includes the 
weak value of a s , {a s ) w = [Y i + Y f +i(Y l x r / )]/(l + r. i - 
Yf) and is derived by the straightforward extension 
of ref. [20] ■ We obtain the six linear equations by 
inserting the controlled parameters in Step 6. Then, 
we can decide the six unknown parameters [guv]- 

The six linear equations can be expressed by £ = A^ 1 ^, 
where 



S • \9xxi 9yyj 9zz} 9xyi Qxzi 9yz) i 

1 + Y iik ■ Yf )k 



C := {Ex(q fe • CTp) - qfc • Pfe}- 



2St k 



[A] j:k := (p fe x q k )u(Yi,k +r/,fc) y 

+ {qfc - p/c(qfe ■ Pfc)]>(r ii/£ x Tf ik )v 



(5) 
(6) 

(7) 



Pt 
Pp 



Hint 


1 Pt 


R(9) 


- 


> 


1 


l_ 








R(q) 


~ 


> 





rotation measurement 



Fig. 1: Quantum circuit representation of our proposal. R(8) 
and -R(q) express the rotation operation of the spin for an arbi- 
trary angle 8 and the measurement direction q for some direc- 
tions, respectively. The dashed-line area indicates the proce- 
dures of the quantum state tomography by changing the angle 



Here, the suffix k denotes the subscript labeling on the 
experimental data and T denotes the matrix transpose. 
Q := {(jl,v)} = {{x,x),(y,y),(z,z),(x,y),(x,z),(y,z)} 
corresponds to the suffix {j} = {1,2,3,4,5,6} since [g^] 
is the symmetric tensor. Therefore, A is the square matrix 
of order six. Our proposal is summarized in Fig. [T] when 
the single spin Hamiltonians, H t and H p , are not con- 
sidered because of Step 5. Compared to the established 
schemes of the quantum process tomography, this advan- 
tage is to evaluate the specific spin-spin interaction with- 
out knowing the relative position between spins |21j , using 
the long-time dynamics [22], and using the Bell measure- 
ment |23| . Furthermore, we would like to emphasize that 
the operations to the target and probe systems are differ- 
ent due to the linearity of the weak measurement while we 
have to use the same operations to the target and probe 
systems from the viewpoint of the conventional quantum 
process tomography. 

Example in Nitrogen Vacancy Center in Dia- 
mond. — In the final part of this Letter, as an example 
to show the reliability of our proposed scheme, we apply 
this to evaluate the hyperfine interaction in a system of a 
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From the example of Table [TJ we obtain the evaluated 
Hamiltonian, 



le-04 le-02 le+00 le-04 le-02 le+00 

8 1 [|xs] 8 1 [fls] 



Fig. 2: Dynamics of the probe spin. Under the controlled 
parameters of Table [l] in which (a)-(f) correspond sequen- 
tially from the top, the dynamics of the expectation values 
of the probe spin is numerically shown. The vertical axis is 
AEx(q -Up), which is defined as the absolute value of the dif- 
ferences between the right-hand side of Eqs. (J3] |4]), i.e., the 
second- and higher-order correction, and the horizontal axis is 
the time scale St of the order of micro seconds. The time at 
which AEx(q-ap) is small has to be chosen at marked positions. 
Since we have not analytically obtained the condition to appear 
in the dent of AEx(q ■ a p ), we now search the points by hand. 
In practical, these points have to be picked up from many tri- 
als changing the set of the controlled parameters. Since the 
interaction Hamiltonian is uniquely determined, we can find 
out these by the above process. 



NV center in diamond [2T] under the above assumption q 
Here, the target and probe spins are the NV electron and 
nearby 13 C nuclear spins, respectively. Because of Step 
5, the single spin Hamiltonians can be ignored; we as- 
sume H t = H p = in the numerics, i.e., q = q. Under 
a given interaction Hamiltonian, we numerically calculate 
the time evolution from the given initial state. Then, we 
calculate the final state for the target spin and the ex- 
pectation value of the probe spin. By Eqs. © - (J7]), 
we numerically evaluate the interaction Hamiltonian. We 
compare the given and evaluated interaction Hamiltonians 
by the error defined as the mean value and the standard 
deviation between the given Hamiltonian [g^ v ] and the 
evaluated one [g^v]'- error := (g) ± a g , where 



id) := Yl (&«'- VO/ 6 ) 



(8) 



E {(^-<^)-<5>}75. (9) 



((i,i/)e(! 



4.98 -6.29 -2.92 
-6.29 4.21 -2.30 
-2.92 -2.30 8.35 



(10) 



2 Since a position of the 13 C nuclear spin cannot be directly ob- 
served in current technology, our proposal seems to be essentially 
needed. 



It is noted that how to decide the choice of the esti- 
mated parameters is not mentioned in this Letter. How- 
ever, they can be in principle picked up from many trials. 
Our proposed scheme can reproduce the original hyper- 
fine interaction with high precision since the evaluated er- 
ror, 0.022 ± 0.063 MHz, is within 1%. It is noted that 
the parameters, ri,p, and q, are determined from Fig. [2] 
Therefore, our proposal may be experimentally realizable 
to evaluate the spin-spin interaction. 

Concluding Remark and Outlook. In conclu- 

sion, we have proposed the procedures from Step 1 to Step 
7 to experimentally evaluate the specific spin-spin inter- 
action ^) in a considered two-qubit system from eq. ((3]) 
obtained in the weak measurement technique. This gives 
the practical advantage to the weak measurement / weak 
value analysis. The weak measurement can be also taken 
as an example of the quantum process tomography and 
quantum estimation. 

There are remaining problems on our proposal. First, 
we have not adjusted our proposal under realistic exper- 
imental constraints. As the first step, we roughly con- 
sidered the realistic experimental setup in the NV center 
in diamond [24]. Second, we have not yet given an opti- 
mized method to choose the controllable parameters while 
the abstract estimation limit is shown |25| . which is re- 
lated to the channel-parameter estimation problem [25] . 
Third, we have not mentioned the scalability [3] of our 
proposal. While it may be possible to overcome the scala- 
bility under the assumption of the spin configuration like 
ref. [5], there are many varieties on this problem by the 
considerable physical qubit. Finally, we did not consider 
the decoherence effects during this scheme. According to 
ref. [12] , it can be included in the decoherence effect for the 
weak value. However, this can be ignored since our pro- 
posal only uses the short-time interaction. This might be 
a practical advantage of our proposal compared with the 
established schemes. We believe that weak measurement 
technique will be a powerful tool to analyze the fundamen- 
tal properties in condensed matter systems and quantum 
optics like our study, despite its naming. 
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